Spin-Blockade in Single and Double Quantum Dots in Magnetic Fields: 

a Correlation Effect 
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The total spin of correlated electrons in a quantum dot changes with magnetic field and this effect 
is generally linked to the change in the total angular momentum from one magic number to another, 
which can be understood in terms of an 'electron molecule' picture for strong fields. Here we propose 
to exploit this fact to realize a spin blockade, i.e., electrons are prohibited to tunnel at specific values 
of the magnetic field. The spin-blockade regions have been obtained by calculating both the ground 
and excited states. In double dots the spin-blockade condition is found to be less stringent than in 
single dots. 
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The Coulomb blockade is one of the highlights in the 
transport properties of mesoscopic systems such as quan- 
tum dots. This is a combined effect of the discreteness of 
energy levels and the electron-electron interaction (charg- 
ing energy). Now, it has recently been suggested that, 
if the total spins of the ground state of N and (N — 1)- 
electrons differ by more than 1/2, the dot is blocked with 
the corresponding peak in the conductance missing at 
zero temperature. This is called the spin-blockade 
and has been studied theoretically for weak electron in- 
teraction regimes. There the Hund's coupling picture, in 
which electrons are accommodated in one-electron states 
with high spins for degenerate states, tends not to realize 
the spin-blockade condition, so that some modifications 
such as an anharmonicity in the confinement potential |^] 
have to be introduced. 

When quantum dots are placed in strong magnetic 
fields, the ground states are known to change dramati- 
cally into the magic- number states |^,^ . This comes from 
the electron correlation effect, since the magic numbers 
for the total angular momentum arise from a combined 
effect of the electron correlation and Pauli's principle, 
persisting even when the Zeeman energy is completely 
ignored. The total angular momentum of the ground 
state jumps from one magic number to another as the 
magnetic field B is varied. 

An important hint that electron correlation is really 
at work is the fact that the total spin (5), where = 
S{S + 1), of the ground state, which dominates how the 
electrons correlate, changes wildly as shown in Fig. |^. 
This happens when the typical Coulomb energy is much 
greater than the single-electron level spacing, where elec- 
tron molecule are formed. In this sense this is genuinely 
an electron-correlation effect — electron correlation has 
been known to dominate the spin states in ordinary cor- 
related electron systems such as the Hubbard model, but 
the present case is a peculiar manifestation in strong 
magnetic fields. 

In the present paper we propose to utilize this elec- 



tron correlation effect to realize a spin blockade. We 
have numerically studied the ground and excited states 
of single dots that contain three or four electrons with 
a parabolic confinement potential and find that the spin 
blockade should indeed be observed. Physically, a key 
observation starts from the fact that the correlated elec- 
tron states in the dot may be thought of as 'electron 
molecules' ||^, which in turn enables us to interpret 
the spin wavefunctions taking part in the spin blockade 
as spin configurations in molecules, which include the 
resonating valence bond (RVB) states, that are usually 
invoked for lattice fermions. We further show that the 
spin-blockade condition is easier to satisfy in double dots 
which can be tuned by controlling the layer separation 
and the strength of the interlayer tunneling. 

So let us start with looking at the total angular mo- 
mentum (L) of 2D electrons confined in a quantum dot 
in a magnetic field , which has a one-to-one correspon- 
dence with the spatial extent (cx of the wavefunc- 
tion. Thus the presence of magic L values signifies that 
the total Coulombic energy of the interacting electrons, 
although roughly a decreasing function of L as the elec- 
trons move further apart for larger L, is not a smooth 
function of the size of the wavefunction, so that jumps in 
L are accompanied by jumps in the size of the wave func- 
tion For example, the total angular momentum L of 
three spin-polarized electrons changes 3 ^ 6 — > 9 — > . . . 
with increasing magnetic field. 

Recently one of the authors has explained this as an 
effect of correlation in the electron configuration, where 
Pauli's exclusion principle dictates group-theoretically 
the manner in which the quantum numbers should ap- 
pear 1^. There, the picture of the 'electron molecule', 
in which the electrons with a specific configuration (tri- 
angle for three electrons, square for four, etc) rotating 
as a whole has turned out to be surprisingly accurate. 
This continues to be the case for larger numbers of elec- 
trons 0]. 

When one considers the spin degrees of freedom, the 
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magic L values are linked with the total spin. This is 
already apparent in the first numerical study of spin 
dependent correlation in quantum dots These 
molecules are characterized by a quantum number, fc^, 
where the spin wave function 5'spin is transformed to 
exp(— 27rfcsi/m)^'spin under the rotation of 27r/m for an 
m— fold symmetric molecule. Then the criterion for the 
magic number, modulo m, reads L + kg = 0(m/2) for m 
odd (even). 

To actually obtain the spin states numerically for dif- 
ferent numbers of electrons, let us consider single GaAs 
quantum dots with three or four electrons in a parabolic 
potential. The electron motion is assumed to be com- 
pletely two dimensional. The Hamiltonian for a single 
dot is 7^ = + Wc, where 
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is the single-electron part, while 
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represents the Coulomb interaction. Here the Hamil- 
tonian is written in second quantized form in a Fock- 
Darwin ^ basis, and e„£^ = (2ri + I + \i\)n {ijjI/4: + 
ojq)^/^— £?ia;c/2— 5*/i£,i?Sz- The dielectric constant is e, 
fiLijQ represents the strength of the parabolic confinement 
potential, Uc = eBlm*c is the cyclotron frequency, m* 
is the effective mass, /i^^ is the Bohr magneton, g* is the 
effective g-factor and is z-component of the spin of a 
single electron. 

We use the confinement potential Tiloq = 6.0 meV. This 
is a little larger than usually estimated values (2 ~ 4 
meV) and is deliberately chosen to reproduce the addi- 
tion energy spectrum The fact that calculations 
with a 1/r interaction require a larger confinement en- 
ergy to reproduce experimental results is considered to 
be a consequence of the modification of the interaction 
potential in real dots |14|. 

In our numerical calculations we have used enough 
states (including higher Landau levels) in the basis to en- 
sure convergence of the ground-state energy within 0.1%. 
Three lowest excited states are also calculated for each 
value of B, which turn out to be the lowest-energy states 
having different angular momenta in the present case. 
Excited states are also obtained with a typical accuracy 
of < 0.1% for an TV = 3 single dot at B = 5 T. 

The total angular momentum, and the total spin of 
the ground state for three- and four-electron systems 
plotted in Fig. |l|, we can see how the magic L val- 
ues go hand in hand with S{N) for N electrons, where 
— S{S + 1) while the z component of S is aligned 



to B: As the magnetic field increase, the ground state 
changes as {L,S) = (1,1/2) ^ (2,1/2) ^ (3,3/2) for 
TV = 3, [L,S) = (0,1) ^ (2,0) ^ (3,1) (4,0) ^ 
(5,1)-^ (6,2) for Ar = 4. 

If we then plot the difference in the total spin, 5(4) — 
5(3), against the magnetic field in the bottom panel of 
Fig. 1^, the spin blockade condition. 



|5(iV)-5(iV-l)|>i 



(3) 



is indeed fulfilled: S jumps from 3/2 to in the region 
4.96 <B < 5.18 T. 

From the magic-number criterion the state with [L = 
2 + 4 X integer) has to have the quantum number ks = 0, 
while the (L = 4 x integer) state kg — 2 for = 4. We 
can make an intriguing identification, by looking at the 
spin density correlation function, that (L, S) — (2, 0) is 
an RVB" state while (4, 0) is an RVB+, where the RVB's 
are defined, for a four-site cluster, as 
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with 
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14-/ I it)) being the spin-singlet 
pair in the electron molecule. The difference of RVB+ 
from RVB~ is that the former lacks the Neel components 
(the last two terms in RVB~) and has the extra phase 
factor -1 for 7r/2 rotation. Although what we have here 
is totally different from lattice fermion systems such as 
the Hubbard model for which RVB is usually conceived, 
the electron-molecule formation has brought about such 
spin configurations. 
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FIG. 1. The total angular momentum and total spin of 
the ground state for a single quantum dot with A'^ = 3 (left) 
or A'^ = 4 electrons (right). The confinement potential is 
assumed to be parabolic with huio ~ 6.0 meV. 
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In this region, the conduction, blocked at zero temper- 
ature, has to occur through an 5 = 1 excited state for 

= 4 at finite temperatures. If the excited states are 
well separated in energy (> 0.1 meV, typical experimen- 
tal resolution for Coulomb diamonds) from the ground 
state for both of the {N — 1)— electron and N— elec- 
tron states, the spin blockade should be observed in the 
Coulomb diamond, which is the differential conductance 
plotted in the plane of source-drain voltage and gate volt- 
age. We have calculated the three lowest excitation en- 
ergies and their total spins for the N = 4 quantum dot 
in Fig. ||. The lowest excited states for iV = 3 and for 
iV = 4 both lie about 0.06meV above the ground state 
around i? = 5.1 T in the spin-blockade region. We can 
make this separation larger (^0.1 meV) for stronger con- 
finement potentials (e.g., 0.09meV around B = 7.4 T for 
Hluq — 8.0 nieV). Such confinement potentials may be 
realized in a gated vertical quantum dot [ p^ . 

The link between the magic L and total S and subse- 
quent spin blockade appears for other numbers of elec- 
trons as well, e.g., between (L, S) = (2, 0) state for iV = 2 
and (6,3/2) state for = 3 for 14.1 < B < 14.8 T. 
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FIG. 2. Top; Excitation energies for A*' = 4 single dot. 
Middle: The same for = 3. Bottom: Difference in the total 
spin, 5(4) — S{3), between the ground state for N = 3,4). 

Now we move on to the double dots, where dots are 
separated in the vertical direction with their centers 
aligned on a common axis. We assume the same con- 
finement potential for the two dots for simplicity. Here 
electrons are Coulomb-correlated both within each layer 
and across the two layers, in the presence of the inter- 
layer tunneling. Recent advances in semiconductor fab- 
rication techniques have enabled fabrication of double 
dots in vertical, triple-barrier structures on submicron 



scales ll^. Theory for the double quantum dots has 
been developed |16|-|2l||, where intriguing features such 
as magic-number states intrinsic to double dots, or a 
singlet-to-triplct spin transition for two-electron system 
have been found ||l^-[T^. 

The Hamiltonian now contains the tunneling term. 



^SAS (ci^^Cnl- + cJ^f_C„£+) , (6) 



while the Coulomb part is now the matrix element of 
e^/e|ri — T2\ for intra-layer interaction, and e^/e(|ri — 
+ (i^)^^^ for inter-layer interaction. The basis is 
\niiiaiai) , where a = ± is an index specifying the two 
dots. 

Thus a double dot is characterized by the parabolic 
confinement potential layer, Tiwo, the layer separation, rf, 
and the strength of the inter-layer tunneling (measured 
by AsaSj the energy gap between the symmetric and an- 
tisymmetric one-electron states). Here we have adopted 
realistic values of Tiujq = 6.0 meV, 10 < d < 50 nm, 0.2 
< AsAS < 2.0 meV. We can now plot in Fig. ^ how high- 
spin states appear on the AgAS — d plane. A high-spin 
state is indeed seen to appear in the upper left region 
of each panel for B = 5.0 T. In the shaded region of 
the right panel, the difference between the total spins is 
S{3) — S{2) = 3/2, fulfilling the spin blockade condition. 
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FIG. 3. The ground-state spin of the double dot for 
A^ = 2(left) or N = 3(right) for B = 5 T for a parabolic 
confinement potential hcjo = 6.0 meV. Shaded region corre- 
sponds to a transition from S = ^ 3/2 

We now focus on a typical point in the shaded region, 
d = 16.0 nm and Asas = 1-2 meV. In Fig. | the total 
energy, total angular momentum and the total spin of 
the ground state for A^ = 2, 3 are plotted. 

The difference between the total spin of two and three 
electrons systems is shown in the bottom panel of Fig. ||. 
The spin-blockade condition is satisfied for 4.0 < B < 
9.3T, which is wider than for the single dot. In the bulk 
bilayer fractional quantum Hall (QH) systems, a phase 
diagram on the Asas — d plane has been considered. If 
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we translate the quantities for the dots, we are work- 
ing in the 'two-component' (correlation-dominated) re- 
gion around the QH-non QH boundary in the language 
for the bilayer QH system. This might have some rele- 
vance to the behavior of the double dots. In Fig. ^, the 
excitation energies for the S ~ 1/2,3/2 states are also 
plotted. The excitation energies for both N = 2 and 

= 3 systems are about 0.12 meV (exceed 0.1 meV) at 
B ~ 6.4 T, which is large enough for the spin blockade to 
be observed. We also notice a level crossing between the 
second and the third excited states around B = 6.9 T for 

= 3, which should appear in the Coulomb diamond. 
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FIG. 4. The total angular momentum and total spin of the 
ground state for double dots with N = 2(left) — 3(right) 
electrons, hujo = 6.0 meV. d = 16.0 nm, and Asas = 1-2 
meV. 
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FIG. 5. Top (middle): Excitation energies for A'^ = 3 
{N = 2) double dots. Bottom: The difference, S(3) - S{2), 
in the total spin for N = 2 and A' = 3 double dots. 



In summary, we have shown that in both single and 
double dots, a spin blockade should occur in some mag- 
netic field region, as an effect of the total spin dominated 
by the magic angular momenta. We would like to thank 
to Seigo Tarucha and Guy Austing for a number of valu- 
able discussions. 
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